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CLAIMS 



I claim: 



1 . A method for imaging prestack seismic data, comprising: 

calculating an individual reflectivity for each frequency in the seismic data; 

calculating a mean reflectivity over the individual reflectivities; 

calculating a variance for the individual reflectivities; 

calculating a variance for the upgoing wavefield, using the mean reflectivity; 

calculating a spatially varying pre-whitening factor, using the variance for the reflectivities 

and the variance for the upgoing wavefield; and 

calculating a reflectivity using the spatially varying pre-whitening factor. 

2. The method of claim 1, wherein the step of calculating a variance for the upgoing wavefield 
comprises applying the following equation: 



of frequencies coj, U(x,ct)j) is the upgoing wavefield, <R(x)> is the mean reflectivity, and D(x,cdj) 
is the downgoing wavefield. 

3. The method of claim 1, wherein the step of calculating a spatially varying pre-whitening 
factor comprises applying the following equation: 



where e(x) is the spatially varying pre-whitening factor, ctu 2 is the variance for the upgoing 
wavefield, and ctr (x) is the variance for the reflectivities. 




where ou 2 (x) is the variance for the upgoing wavefield, x is the spatial location, n is the number 
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4. The method of claim 1, wherein the step of calculating a reflectivity using the spatially 
varying pre-whitening factor comprises applying the following equation: 

|Z)(x,a>,)| 2 + ^ 

where R(x) is the reflectivity, x is the spatial location, n is the number of frequencies coj, U(x,o)j) 
is the upgoing wavefield, D(x,coj) is the downgoing wavefield, ou 2 is the variance for the 
upgoing wavefield, and a R 2 (x) is the variance for the reflectivities. 

5. A method for imaging prestack seismic data, wherein a least squares approach comprises 
applying the following equation: 

±£/>*(x,<»,)C/(x,c> y ) 

*(*)=T^ • 

-jTz>*(x,«»,) D(x,o)j) + e 
n j=l 

where R(x) is the reflectivity, x is the spatial location, n is the number of frequencies coj, U(x,ojj) 
is the upgoing wavefield, D(x,a>j) is the downgoing wavefield, and s is a pre-whitening factor. 
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